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We show that all current formalisms for quarks in lattice QCD are consistent in the quenched continuum 
limit, as they should be. We improve on previous extrapolations to this limit, and the understanding of lattice 
systematic errors there, by using a constrained fit including both leading and sub-leading dependence on a. 



1. INTRODUCTION 

Our increasing ability to carry out unquenched 
simulations means that the quenched approxima- 
tion is no longer required. There is an outstand- 
ing issue associated with the quenched continuum 
limit, however, which we address here. At the 
same time we point out the ingredients that are 
necessary to obtain an accurate result in the con- 
tinuum limit, since these ingredients will also be 
important for unquenched simulations. 

The outstanding issue is the question of 
whether all quark formulations give the same re- 
sult in the quenched continuum limit. It is rel- 
atively trivial to demonstrate that they should, 
because this is the limit in which it is easy to 
analyse the formulations. At Lattice 2000, how- 
ever, it seemed not to be the case U|- A plot 
was made of the ratio of nucleon to vector meson 
mass at a fixed physical quark mass, as a func- 
tion of the lattice spacing. The physical quark 
mass was chosen as the point at which the pseu- 
doscalar meson mass divided by the vector meson 
mass was 0.5 (in those days quite a low mass to 
reach). The lattice spacing was given in units of 
the vector meson mass. Because of the instabil- 
ity of vector mesons this is not a plot that allows 
for a precision test against experimental results, 
even if it used a quark mass value from the real 
world and was not in the quenched approxima- 
tion. However, it can be used as a comparison 



of different formalisms, provided the results have 
good statistical precision and are on large enough 
volumes. 

The original plot showed a 6a disagreement be- 
tween the results extrapolated to a = from the 
Wilson formulation |2] and the unimproved stag- 
gered formulation 0. The Wilson results were 
extrapolated linearly in a and the staggered re- 
sults quadratically. Clover results were also in- 
cluded but their variation suggests a statistical 
error that means they are not particularly useful 
so we drop those results in this discussion. 

The disagreement, if true, would represent a 
major problem for lattice calculations and call the 
whole approach into question. With many large 
dynamical quark projects well underway and oth- 
ers planned, it is therefore time to readdress this 
issue, following on from In addition, there are 
more quark formulations in widespread use than 
there were in 2000, so the comparison of results 
can be significantly extended. 

2. NEW ANALYSIS AND RESULTS 

Figure ^ shows an updated plot of the Wil- 
son and unimproved staggered results from Lat- 
tice 2000, along with new results from improved 
staggered (asqtad) quarks on improved glue H] 
and the perfect action 1^1. The new results were 
chosen for small statistical errors, being available 
at the quark mass required here (with interpola- 
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Figure 1. mj^ jmy for mpg/mv = 0.5 for a vari- 
ety of quark formalisms in the quenched approxi- 
mation. The curves are for a fit, described in the 
text, with a single continuum limit, marked with 
a filled circle. The filled grey triangle at a=0 
shows the previous continuum limit for Wilson 
quarks obtained with a purely linear fit. 



tion), for more than one value of the lattice spac- 
ing and with a reasonably large physical volume. 
Note that the x axis is now a? instead of a. 

The lines represent a simultaneous fit to all the 
data with a single continuum limit imposed. A 
good is obtained. The fit allows for leading 
and higher order polynomials in a appropriate to 
each formulation (i.e even powers for all except 
Wilson) with a constraint placed upon the coef- 
ficients to avoid losing control of the fit. The ef- 
fectiveness of this Bayesian approach [7] has been 
widely demonstrated for fitting correlators, chi- 
ral extrapolations etc. - all situations in which 
a critical issue is the systematic error in the final 
result which arises from leaving out higher order 
terms in an expansion. In all these cases we have 
a good physical understanding of the expansion 



Figure 2. my in units of ri at tops''i=0.807. 
Results are from different quark formalisms and 
the fit has a joint continuum limit, as described 
in the text. 



and can therefore place constraints on the higher 
order terms. Taking the continuum limit is just 
such a case. We expect the scale of discretisation 
errors to be set by the size of internal momenta 
inside the hadron. This should be roughly a few 
hundred MeV. Here the priors on the coefficients 
in the polynomial in myo were taken to be ±0.5, 
and 5 terms were included in the expansion for 
each formulation (a constant plus 4 appropriate 
powers) . 

The result for the joint continuum limit is 
1.373(10), 2(7 below the previous unimproved 
staggered continuum limit but well above the 
Wilson one. With hindsight this is not surprising 
because a purely linear fit to Wilson results ne- 
glects terms which are not small for these results. 
We would expect the quadratic terms to appear 
as, say, (mya/2)^ and this is 5% for the coarsest 
lattice used. A systematic error has to be added 
to a purely linear extrapolation to take this into 
account. For the improved formalisms this is not 
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Figure 3. ttin in units of ri at mpsri=0.807. 
Results are from different quark formalisms and 
the fit has a joint contintium limit, as described 
in the text. 



so much of an issue because the terms neglected 
by a leading order extrapolation are C'(TOya/2)'*, 
1% for the improved staggered and perfect actions 
at their coarsest spacing. In all cases, however, a 
Bayesian analysis allowing for higher order terms 
is useful in assessing the systematic errors. 

We now turn to other scaling plots, Figures |21 
and 13 for which a lot more data is available. 
These are plots of mjv and my in units of ri, 
vs (a/ri)^ We include results from Wil- 

son [2], clover perfect action |H] and do- 

main wall quarks jlll as well as a variety of stag- 
gered results with improved and unimproved glue. 
Again a good fit with a single continuum limit 
is readily obtained and the curves for this are 
shown. The results are for a fixed quark mass 
given by tupsTq — 1.127 using r^/ri = 1.397 in 
the quenched continuum limit 'S'. 

The authors would welcome additional results 
for any of these graphs. 



3. CONCLUSION 

We show, reassuringly, that all current quark 
formalisms give the same answer for nucleon and 
rho masses, at fixed arbitrary quark mass, in the 
quenched continuum limit. We used an analy- 
sis of discretisation effects which allows the in- 
clusion of leading and higher order terms in the 
lattice spacing. This is particularly necessary for 
unimproved actions, emphasising once again the 
importance of the improvement programme in re- 
ducing systematic errors in lattice calculations. 
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